Exact solutions for Big Bounce in loop quantum cosmology 
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In this paper we study the flat (fc = 0) cosmological FRW model with holonomy corrections of Loop 
Quantum Gravity. The considered universe contains a massless scalar field and the cosmological 
constant A. We find analytical solutions for this model in different configurations and investigate 
its dynamical behaviour in the whole phase space. We show the explicit infiuence of A on the 
qualitative and quantitative character of solutions. Even in the case of positive A the oscillating 
solutions without the initial and final singularity appear as a generic case for some quantisation 
schemes. 



I. INTRODUCTION 

In recent years Loop Quantum Cosmology (LQC) has 
inspired realisation of the cosmological scenario in which 
the initial singularity is replaced by the bounce. In this 
picture, the Universe is initially in the contracting phase, 
reaches the minimal, nonzero volume and, thanks to 
quantum repulsion, evolves toward the expanding phase. 
Such a scenario has been extensively studied with use of 
the numerical methods [H, Ql- However, as it was shown 
for example in exact solutions for bouncing universe 
with dust and cosmological constant can be found. The 
aim of the present paper is to show that analytical solu- 
tions can also be obtained for the bouncing models arising 
from LQC. The main advantage of such exact solutions 
is that they allow for investigations in whole ranges of 
the parameter domains. 

In this paper we consider the flat FRW model with 
a free scalar field and with the cosmological constant. 
Quantum effects are introduced in terms of correction to 
the classical theory. Generally one considers two types 
of of quantum correction: correction from inverse vol- 
ume and holonomy corrections. The leading effect of 
the volume corrections is the appearance of the super- 
inflationary phase. The effect of holonomy corrections, 
on the other hand, is the appearance of a bounce instead 
of singularity. The aim of this paper is to investigate 
analytically these effects in a flat FRW model. That is 
to say, we neglect corrections from inverse volume, these 
effects however, has been extensively studied elsewhere. 
Moreover, these two types of corrections are not equally 
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important in the same regimes. The inverse volume cor- 
rections are mainly important for small values of the scale 
factor, whereas holonomy corrections are mainly impor- 
tant for large values of the Hubble parameter. In other 
words, when the minimal scale factor (during the bounce) 
is large enough, the effects of inverse volume corrections 
can be neglected. 

The flat FRW model in the Loop Quantum Cosmology 
has been flrst investigated in the pioneer works of Bo- 
jowald 5, 6] and later improved in the works of Ashtckar, 
Pawlowski and Singh (JH, 0|- Bojowald's original de- 
scription of the quantum universe in currently explored 
in the number of works and regarded as a parallel line 
of research [9^,^01. In the present paper, we restrict our- 
selves to the flat FRW models arising in the framework 
proposed by Ashtekar and co-workers. Beside the flat 
models this approach has also been applied to the FRW 
k = ±l models in [1, [O, [H, Hi and Bianchi I in [i3,[l|]. 
In these models the unambiguity in the choice of the ele- 
mentary area for the holonomy corrections appear. In the 
present paper we consider two kind of approaches to this 
problem: the so called /i— scheme and /xq— scheme (for a 
more detailed description see Appendix \Q . We find an- 
alytical solutions for the considered models in these two 
schemes. 

The Hamiltonian of the considered model is given by 



-ffcff = - 



87rG72 



sin (Ac) ] 1 I ,3/2 A 

(1) 

In Appendix [X] we show the derivation of this Hamilto- 
nian in the Loop Quantum Gravity setting. The canon- 
ical variables for the gravitational field are {c,p) and for 
the scalar field {4>,p^). The canonical variables for the 
gravitational field can be expressed in terms of the stan- 



dard FRW variables (c, |p|) = {jaVg^^ ,a'^V^^^'^). Where 
the factor 7 is called Barbero-Immirzi parameter and is a 
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constant of the theory, and Vq is the volume of the fidu- 
cial cell. The volume Vq is just a scaling factor and can 
be chosen arbitrarily in the domain Vq € IR.+ . Since p is 
the more natural variable than a here, we present mostly 
p(t) in the figures, a is always the positive square root 
of p so the shape of the graphs would be essentially the 
same when drawn with a. The equations of motions can 
now be derived with the use of the Hamilton equation 



f = {f,H,s} 
where the Poisson bracket is defined as follows 
87rG7 



(2) 



{f,9} = 



d£dg_ _ d£dg_ 
dc dp dp dc 
df dg 
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df dg 

d(j) dp^ dp^ d(t> 



(3) 



From this definition we can immediately retrieve the el- 
ementary brackets 

{c,p} = ^^ and = (4) 

With use of the Hamiltonian ([T]) and equation ^ we 
obtain equations of motion for the canonical variables 



P 



1 M 



sin (/ic) cos (/ic) , 



sin (/ic) 



sgn(p)^|-^ + sgn(p)^yH: 



'■^ = 0, 



(5) 



where k = SttG. The Hamiltonian constraint iJeff — 
implies 



sin (/ic) 



l\E±. A 
^+ 3- 



3 2br 



(6) 



The variable p, corresponds to the dimensionless length 
of the edge of the elementary loop and can be written in 
the general form 



Mb) = 



(7) 



where —1/2 < n < and ^ is a constant ^ > (this 
comes from the fact that p is positively defined). The 
choice of n and ^ depends on the particular scheme in 
the holonomy corrections. In particular, boundary values 
correspond to the cases when p, is the physical distance 
{n = —1/2, /i— scheme) and coordinate distance (n = 0, 
/LtQ- scheme). However, the n = case does not lead to 
the correct classical limit. When n > 0, the classical limit 



can not be recovered either. Only for negative values of 
n is the classical limit p — > cxd correctly recovered 



lim 



sin {n{p)c) 

Kp) 



(8) 



Strict motivation of the domain of the parameter n 
comes form the investigation of the lattice states f25| . 
The number of the lattice blocks is expressed as, N = 
Vq/Iq where Iq is the average length of the lattice edge. 
This value is connected to the earlier introduced length 
/i, namely M — p,^^{p). During evolution the increase of 
the total volume is due to the increase of the spin labels 
on the graph edges or due to the increase of the num- 
ber of vortices. In this former case the number of the 
lattice blocks is constant during evolution, J\f = const. 
Otherwise, when the spin labels do not change, the num- 
ber of vortices scales with the volume. A/" oc bP^^. The 
physical evolution correspond to something in the mid- 
dle, it means the power index lies in the range [0,3/2]. 
Applying the definition of Af we see that the considered 
boundary values translate to the domain of n introduced 
earlier, n e [—1/2,0]. More detailed investigation of the 
considered ambiguities can be found in the papers p^flGj 
and in the Appendix \X[ 

Combining equations ([6]), ([7]) and the first one from 
the set of equations ([5]) we obtain 
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n,\pr'+^n\p\'~^ni\p\''''' 



where new parameters are defined as follow 



|2n+3 



(9) 

(10) 

(11) 
(12) 
(13) 
(14) 

Equation ([5]) is, in fact, a modified Friedmann equation 

(15) 
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Qui = 


l^'l^pt, 
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H = pH , 

3 3 ' 



where the effective constants are expressed as follow 



Goff — G 
Acff = A 



1-^ 



Pc 

Pc Kpc 



P A 
1-2— 



and 



2br 



(16) 
(17) 

(18) 
(19) 
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We will study the solutions of the equation ([9]) for both 
//—scheme and ^o~schenie. 

The organisation of the text is the following. In sec- 
tion [IT] we consider models with A = 0. We find solutions 
of the equations ([9]) for both /i— scheme and /io~scheme. 
Next, in section [TTTl we add to our considerations a non- 
vanishing cosmological constant A. We carry out an anal- 
ysis similar to the case without lambda. We find analyt- 
ical solutions of the equation([9]) for /I— scheme. Then, 
we study the behaviour of this case in fj,Q— scheme. In 
section fVl we summarise the results. 



II. MODELS WITH A = 

In this section we begin our considerations with the 
model without A. Equation ^ is then simplified to the 
form 



(20) 



We solve this equation for both /x— scheme and 
/LiQ— scheme in the present section. 



The main property of this solutions is that p{t) never 
reaches zero value for non vanishing p^. The minimal 
value of p{t) is given by 



bmini = "fl 



1/3 



(26) 



3 V^pi, 

We show the solution (|25p in the Fig. [TJ The dynamical 



A. /_t— scheme (n : 



-1/2) 



In the /i— scheme, as is explained in the Appendix \^ 
the fi is expressed as 



(21) 



where A = 2^/3TT^lpi is the area gap. So ^ = \/A and 
n = -l/2. 

To solve (l^ni) in the considered scheme we introduce a 
new dependent variables u in the form 




\pit)\=u'/'it) 



(22) 



With use of the variable u the equation ([20)1 takes the 
form 



(23) 



and has a solution in the form of a second order polyno- 
mial 



^ + ^n,t' - Hf^iCit -f- ^Cf, (24) 
sil 4 4 4 



where Ci is a constant of integration. We can choose 
now Ci = 0, so that the minimum of u(t) occurs for 
t — 0. Going back to the canonical variable p we obtain 
a bouncing solution 



p{t) = sgn(p) 



f7i 4 



1/3 



(25) 



FIG. 1: Typical solution with A = in the /i— scheme. The 
canonical variable p is expressed in the [lp{\ units and time t 
in the [Ipi] units. 

behaviour in this model is simple. For negative times we 
have a contracting pre-Big Bang Universe. For t > we 
have a Big Bang evolution from minimal Pmin (|26p . It is, 
however, convenient to call this kind of stage Big Bounce 
rather than Big Bang because of initial singularity avoid- 



B. /^o— Scheme (n = 0) 

In the /iQ— Scheme the jl is expressed as 

3^/3 



M = Mo 



(27) 



So (, — Ho and n = 0. To solve (|20p in this scheme, we 
change the time 



dt = y^\p{u)\du (28) 

and introduce a new variable y as follows 

y{u) = \p{u)\\ (29) 

Applying this new parametrisation to the equation (I20p 
leads to the equation in the form 



(30) 
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which has a solution 



y{u) ^ ^ + fliu^ - 2fliCiu + QiCl (31) 

We can now choose Ci = so that the minimal value of 
y(u) occurs for u — 0. Now we can go back to the initial 
parameters p and t, then 



p{u) = sgn{p)\j'-^ + Q.IV? 



t{u) 



du' 



fliii 



-,1/4 



Introducing the variable 



III 



(32) 
(33) 

(34) 



we can rewrite integral (|33p to the simplest form 




FIG. 2: Typical solution with A = in the /xo— scheme. The 
canonical variable p is expressed in [lp{\ units and time t in 
the [Ipi] units. 



1/4 



S2j 

Jo 



and the solution of such integral is given by 



Jo 



1 3 3 



, (36) 



where 2F1 is the hypergeometric function, defined as 

{ai)k...{ap)k 



oFq[ai, Up; bi, ...,bq; x]=^ 



(6i),...(6,)fc fc! 



(37) 



where {a)k is the Pochhammer symbol defined as follow 
T{a + k) 



{a)k 



r(a) 



a(a+l)...(a + fc- 1). (38) 



This solution is very similar to the one in /i— scheme. 
However, the time parametrisation is expressed in a more 
complex way. We show this solution in the Fig. [21 In 
this case minimal value of p{t) is expressed as 



(39) 



III. MODELS WITH A / 

In this section, we investigate the general model with 
non vanishing cosmological constant. It will be useful to 
write equation ^ in the form 



ni\pr^ + nu\p\' 



^ ^111, |2ri-2 ^^V , ,2n+l 



(40) 



We see that when we perform the multiplication in this 
equation and define 



IV 



ii"iii 



I"V 



(41) 



we recover equation In this and the next section 

we use equivalently Aq — and a = -^p^ to simplify 
notation. 



A. /i— scheme {n — —1/2) 



In this case equation (|40p can be rewritten to the form 



^) =7i|p|^ + 7iibl ^-7iii|p| " 



where the parameters are defined as 



71 = fill -fly^-A 



711 = ill - iliv = -np^ 

o 1 2 2^2 4 

7III = "III = 9^ ^ ^ P't'- 



A 



2^2 



1-2-7^^ 



(42) 

(43) 
(44) 
(45) 



To solve equation (j42p we re-parametrise the time vari- 
able 

dt = \p{u)\^du, (46) 
and introduce a new variable p as follows 

y{u) = \p{u)\\ (47) 
This change of variables leads to the equation in the form 



^) =W(2/-yi)(y-y2), (48) 
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where 



2/1 



y2 = -7T 



1 ^P^T^C^ 

6 [1 - A^2^2] 

2 A ■ 



(49) 
(50) 



There are three general types of solutions correspond- 
ing to the values of the parameters (71, yi,y2)- We sum- 
marise these possibilities in the table below. 
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< 


> 


> 




> 
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< 


71 


< 


> 


< 


yi 


> 


> 


< 


y2 


> 


< 


< 



It is important to notice that the product T = 71 • yi • 

2 

2/2 — "12^^ is negative in all cases. This property will 
be useful to solve equation In the Fig. [3] we show 

values of the roots yi and 2/2 as functions of A. Thus, 



0.4 



0.1 


-0.1^ 
-0.2: 
-0.3 

-0.4 



y2 


yi 




1 


/ ^ 
y2 


yi 



50 100 
A 



150 



200 



FIG. 3: Parameters yi and 2/2 as functions of A. Region 
1 corresponds to oscillatory solution, region 2 to bouncing 
solution. There is no physical solutions in region 3. The 
parameter A is expressed in units of [mpi]. 



there are two values of cosmological constants where the 
signs of the roots change, namely A = and 



A, 



0,p. 



'PI 



27r7'^ 



133.4mf,i, 



(51) 



where we have used the value of the 7 = ln2/(7r\/3) 
calculated in the work More recent investigation of 
the black hole entropy indicate however that the value 
of the Barbero-Immirzi parameter is ln2/7r < 7 < ln3/7r 
19 1. In particular, Meissner has calculated jm — 0.12738 
20|. However, this freedom in the choice of 7 does not 
change the qualitative results. Only the region of the 



parameter space where the particular kind of motion is 
allowed can be shifted. 

We now perform a change of variables in equation 
in the form 



w{u) 



and 



We also introduce the parameters 
1 1 

Wi = , W2 = , 

yi y2 

and then equation (I48p takes the form 
dw ^ ^ 



— — -(wi — w)(w2 — w). 
dX, 



(52) 
(53) 
(54) 

(55) 



The equation (j55p is the equation of shifted harmonic 
oscillator and its solution is 



w{x) 



6a 



1 - 2 



A 



- COS (x) 



(56) 



where we set the integration constant to zero. It is now 
possible to go back to the initial parameters p and t which 
are expressed as follows 



p{u) = sgn(p) 



6a 



(57) 



t{u) 



1- 24-+ cos (3^.) 
du'p^iu), (58) 



The integral ([55)1 can be easily solved but solutions de- 
pend on the parameters in ((57)) . 
We define now an integral 



I{x, a) 



dx' 



Iq a + cos x' 
so the expression (j58|) can now be written as 



tiu) 



3 a/it 



(59) 



(60) 



Solutions of integral ([55)1 depend on the value of param- 
eter a and read 



/(a;,|a|>l) = 

Va — 1 

I{x,a = 1) = tan (^^^ , 
l{x,0< \a\ < 1) = - 

I{x, a = 0) = In 



: arctan 



(a - l)tan(f) 



1 



larctanh 



(a - l)tan(f) 



— tanx 



cos X 
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-TO _8 -6 -4 -2 2 4 6 8 10 
u 

FIG. 4: Parametric solution with A < 0. Oscillating curve 
(green) is 1^(11)1 and the increasing curve (red) is t{u). The 
canonical variable p is expressed in the [Ipi] units and time t 
in the [Zpi] units. 



-4 -8 -6 -4 -2 2 4 6 8 10 
u 

FIG. 6: Parametric solution with < A < Ao.^. Top curve 
(green) shows \p{u)\ and bottom curve (red) is t(u). The 
canonical variable p is expressed in units of [Ipi] and time t in 
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FIG. 5: Solution |p(f)|with A < 0. The canonical variable p 
is expressed in units of [Ipi] and time t in [Ipi]. 



7 

6^ 
5 

P 4 
3 

2- 
1: 






-10 



10 



15 



FIG. 7: Solution \p{t)\ with < A < Ao,p.. The canonical 
variable p is expressed in units of [Ipi] and time t in [Zpi]. 



In Fig. [4] we show the solution for A < - or equiva- 
lently with parameter a > 1. 

In Fig. Owe show the solution for < A < Ag^^ which 
corresponds to < |a| < 1. 



B. po — Scheme (n — 0) 

In this section we study the last case in which n = 
and ^ = 3-\/3/2. Equation (|40|) can then be written in 
the form 



= \p\-' {n, + nu\p\')w{\p\), (61) 



where 



W{\p\) = \p\\3 + \p\\2 + \p\ai + ao, (62) 



and the polynomial's coefficients are expressed as 



1 



ao = - 

ai = 0, 

0-2 = 1, 

03 ■ 



2 2/-2 

:i^P4.1 ^ 



3a 



A 2^2 

3 ' ^ 



A 



The discriminant of polynomial (j62p is 



(63) 

(64) 
(65) 

(66) 



A 



a^al — 40^ as — 4a2ao + 1800010203 — 27a§a: 



(67) 



Inserting the values of parameters {00,01,02,03} listed 
above we obtain 



A 



A, 



12a / 3a A 

- 27 



A, 



(68) 
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For A > 0, or 



81 



> 



A, 



A 



A, 



(69) 



polynomial ([5^ has three real roots. When relation 
is fulfilled and A > 0, oscillatory solutions occur. For A < 
equation (j6ip has bouncing type solutions. This can be 
seen when we redefine equation to the point particle 
form. This approach is useful for qualitative analysis 
and will be fully used in the next section. However, here 
we will use it to distinguish between different types of 
solutions. 

Introducing a new time variable 



dt = 



\P\ 



3/2 



zdu. 



(70) 



we can rewrite equation (1611) to the form 



n 



1 



V{p) = E 



a 



2Ao,, 



where the potential function 
1 



V{j>) 



A 



2'P''('- A 



(71) 



(72) 



Equation ([71]) has the form of Hamiltonian constraint 
for a point particle in a potential well. We see that for 
A < potential ([7^ has only one extremum (for p = 0) 
and only a bouncing solution is possible. For A < 
potential (|72p has a minimum for jpn 



_ 2 ^^0,1^0 



In 

this case physical solutions correspond to the condition 
E > V{\pniin\) and the energy of the imagined particle in 
the potential well is greater than the minimum of well. 
The particle then oscillates between the boundaries of 
the potential. The condition E > y(|pmin|) is equivalent 
to relation (|69p calculated from the discriminant. 
Upon introducing the parameters 



1 



1 



92 := Y2"2 ~ 4"3ai, 



1 

and the variable v 

4 

bl 



216 



1 2 



03 



Q2 



(73) 
(74) 

(75) 



then equation (l6T|) takes on the form of the Weierstrass 
equation 



92V - 53. 



(76) 



appendix to the article [l7|. The solution of equation 
(ffS)) is the Weierstrass p-function 



v{u) = p{u - uo;g2,53), 



where 



92 



93 



1 

12' 



1 
216 



A 



16 V-'^O./io/ ^0,^0 



3a 



So the parametric solution for the parameter p is 

"1 

3 



p{u) = sgii{p)^^^^ 



A 



4p(u-uo;52,5'3) 



(77) 

(78) 
(79) 

(80) 



The time variable can be expressed as the integral 

/ \ 3/2 

t(u) = / fh,' ^ (81) 



f7i + f7„N|«(u')-i|: 



In the Fig. [8] and [9] we show an exemplary parametric 
bounce solution with A < and a possible oscillatory 
solution with A > 0. 




FIG. 8: Parametric solution \p{u)\ with A < 0. The canonical 
variable p is expressed in units of [ipi], u is dimensionless. 

In general the solution is expressible as an explicit func- 
tion of time, by means of the so called Abelian functions. 
However, we chose not to employ them here as the appro- 
priate formulae are not as clear, and the commonly used 
numerical packages do not allow for their direct plotting. 
The fact of existence of such solutions allows us to as- 
sume that the above integral is well defined, and so is 
the solution itself. 



IV. QUALITATIVE METHODS OF 
DIFFERENTIAL EQUATIONS IN STUDY 
EVOLUTIONAL PATHS 



Detailed analysis and plots of solutions of this equation 
for different values of parameters can be found in the 



The main advantage of using qualitative methods of 
differential equations (dynamical systems methods) is the 
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FIG. 9: Parametric solution \p{u)\ with A > and relation 
(|69l) satisfied. The canonical variable p is expressed in units 
of [Ipi] and u is dimensionless. 



investigation of all solutions for all admissible initial con- 
ditions. We demonstrate that dynamics of the model can 
be reduced to the form of two dimensional autonomous 
dynamical system. In our case the phase space is 2D 
{p,p). First we can find the solutions corresponding to 
vanishing of the right hand sides of the system which 
are called critical points. The information about their 
stability and character is contained in the linearisation 
matrix around a given critical point. In the considered 
case, the dynamical system is of the Newtonian type. For 
such a system, the characteristic equation which deter- 
mines the eigenvalues of the linearisation matrix at the 
critical point is of the form -f- d^V{p)/dp^ where V{p) 
is a potential function and \p\ — o^Vq^^, a being the 
scale factor. As it is well known for dynamical systems 
of the Newtonian type, only two types of critical points 
are admissible. If the diagram of the potential function 
is upper convex then eigenvalues are real and of opposite 
signs, and the corresponding critical point is of the sad- 
dle type. In the opposite case, if d^V{p)/dp^ > then 
the eigenvalues are purely imaginary and conjugate. The 
corresponding critical point is of the centre type. 
Equation can be written in the form 

(82) 

where a = ^p'^ or equivalently as 

1 l?l|-2"+3 A 



4a + A^^ =P =1^1 



-i'e'-^\p\)-i'ea, 

(83) 

where we have made the following time reparametrisation 
t — s- M 



d_ _ 1 \p\'''+^/^ d 
du~2 /^^A^p^sdf 



(84) 



Now we are able to write the Hamiltonian constraint 
in the form analogous to the particle of the unit mass 



moving in the one dimensional potential well 



1 



Vip) = E 
where the potential function 



2<-2 ' 



2' 



Vip) 



1 



-2n 



2/-2 



A, 



(85) 



(86) 



As we can see, the constant — 7^^^q;/2 plays the role of 
the total energy of the fictitious particle. The domain 
admissible for motion in the configuration space is deter- 
mined by the condition V{p) -\- 7^^^ck/2 < 0. 

A dynamical system of the Hamiltonian type has the 
following form 



P =y. 

v' = 



(87) 



where the prime denotes differentiation with respect to a 
new re-parametrised time variable which is a monotonous 
function of the original, cosmological time. 

The structure of the phase plane is organised by the 
number and location of critical points. In our case, crit- 
ical points in the finite domain are located only on the 
line p' = y = and the second coordinate is determined 
form the equation —dV/dp — which is 



A, 



0. 



The number of critical points depends on the value of A 
and n. We can distinguish two cases: 



• for -1/2 < 71 < 0: 
A<0: i\p\,p') 
A>0: (\pIp') 

i\plp') = 



(0,0); 
(0, 0) and 

2„ + l/ 2(l-Tl) " 



• for n 



-1/2: 

3 



A^;^: (b|,p') = (0,0); 



A: 



degenerate V{p) 



0. 



The full analysis of the behaviour of trajectories re- 
quires investigation also at the infinity. To this aim we 
introduce radial coordinates on the phase plane for com- 
pactification of the plane by adjoining the circle at infin- 
ity P = TZ7 cos e, y ^ -^sinO. 

The phase portraits for both cases are shown at Figs. 
[To] and [m 

For the general case of A > we can write the para- 
metric equation of the boundary a = of the physically 
admissible region in the phase space, namely 

p2 = lA|pp(l-7^e^|br+i), (89) 

where the dot denotes differentiation with respect to cos- 
mological time t. This equation greatly simplifies for the 
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special case n = —1/2, and we receive the value of the 
Hubble parameter at the boundary 



(90) 



In this special case a = and rt = we can integrate 
eq.(l) for A > and choosing the integration constant 
equal to zero 



p{t) = sgn(p) 



1 



^ cosh2(yAt)' 



(91) 



which gives the maximal value o parameter p 

Ao 



\Pn 



A 



(92) 



In the same case with A < and choosing integration 
constant equal to zero we have 



p{t) = sgn{p) 



A, 



1 



A 



sinh 



(93) 



For a — and A < Ao,/i in the case n = —1/2 we 
obtain the de Sitter solution 



Pit) =. sgn(p) const . exp |2y 1 - ^ ^^tj . (94) 

These solutions represent the lines on the boundaries 
of the physically admissible domains. 



V. SUMMARY 

We have studied dynamics and analytical solutions 
of the flat Friedmann-Robertson- Walker cosmological 
model with a free scalar field and the cosmological con- 
stant, modified by the holonomy corrections of Loop 
Quantum Gravity. 

We performed calculations in two setups called 
/2— scheme and /.jq- scheme, explained in the appendix A. 
We have explored whole A G R range and whole allowed 
p| e M+ U {0} range. In the case of ^—scheme resulting 
solutions are of oscillating type for A G ]R_ and bounc- 
ing type for A G (0,Ao). For A G [Ao,oo) there are no 
physical solutions. In the case of /xq- scheme for A G M_ 
bouncing solutions occur. When both A G K.+ and rela- 
tion (|69p are fulfilled, oscillatory behaviour occurs. Oth- 
erwise, bouncing solutions appear. In all considered cases 
with p^ G M+ the initial singularity is avoided. 

We have investigated the evolutional paths of the 
model, from the point of view of qualitative methods of 
dynamical systems of differential equations. We found 
that in the special case of n = — 1/2 the boundary trajec- 
tory (a = 0) approaches the de Sitter state; and demon- 
strate that in the case of positive cosmological constant 



there are two types of dynamical behaviours in the finite 
domain. For the case of n = there appear oscillating 
solutions without the initial and final singularities, and 
that they change into bouncing for > n > —1/2. 

The results of this paper can give helpful background 
dynamics to study variety of physical phenomenas during 
the bounce epoch in Loop Quantum Cosmology. For ex- 
ample, the interesting question of the fluctuations like 
gravitational waves f2li l23j or scalar perturbations 
[24| during this period. We tried to show that numer- 
ical calculations can be "shifted" one step further, since 
the basic model is explicitly solvable, and can be treated 
as starting ground for more complex problems, like the 
above, which cannot be solved analytically. 
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APPENDIX A: FLAT FRW MODEL AND 
HOLONOMY CORRECTIONS IN LOOP 
QUANTUM GRAVITY 

In this appendix, we derive the form of the Hamilto- 
nian ([1]) considered in the paper. 

The FRW fc = spacetime metric can be written as 



(Al) 



where N(x) is the lapse function and the spatial part of 
the metric is expressed as 



(A2) 



In this expression "g^^ is fiducial metric and are co- 
triads dual to the triads °e", °w*(°ej) = (Sj where °a;* = 
"w^dx" and °ei = °e^da- From these triads we construct 
the Ashtekar variables 

Al Tl+^Kl^c"u;l, (A3) 

Ef EE ^/\d^\e^ ^pV^'^el (A4) 

where 

\p\ = a\ (A5) 

c = 7d. (A6) 

Note that the Gaussian constraint implies that p ^ ~p 
leads to the same physical results. The factor 7 is called 
Barbero-Immirzi parameter. In the definition (jA3P the 
spin connection is defined as 



ri 



1 



(A7) 
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-0.7 5 



V[p] 



V[p] 



0.5 0.5 



1.5 



FIG. 10: The phase space diagram for the case n = and a) A > and b) A < 0. The physical domain admissible for motion is 
shaded. Note that for the case (a) the boundary of admissible for motion is bounded by a homoclinic orbit and all solutions in 
this area are oscillating without initial and final singularities. All trajectories situated in physical region posses the minimum 
value of scale factor. 



and the extrinsic curvature is defined as 
1 



2N 



(A8) 



which corresponds to iiT* := Kabe-\. 

The scalar constraint, in Ashtekar variables, has the 
form 



1 



(A9) 



167rG7s \/\AetE\ 
- 2(l + 7')iffX, > 
where field strength is expressed as 

F^, = daAl ~ dbAl + e%J^,A', . (AlO) 



With use of and (|AlOl) the Hamiltonian 

assumes the form 



He 



3^0 



'\P\i? 



(All) 



where we have assumed a gauge of N{x) = 1. The con- 
stant Vb is the volume of the fiducial cell. This volume 
can be chosen arbitrarily. It is convenient to absorb the 
factor Vb by redefinition 

p^pV^'\ c^cy^'\ (A12) 
The holonomy along a curve a e S is defined as follows 



ha = "Pexp / TiA\dx°-, 



(A13) 



where 2iTi = ai and Ui are the Pauli matrices. From 
this definition we can calculate holonomy in the direction 

°e°;da and the length ^JLVQ^^ 



hi = exp 



1/3 







icos I y) + 2r^sin (^y 



(A14) 
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V[p] 



1.5 



V[p] 



1.5 



FIG. 11: The phase space diagram for the case n = —1/2 and a) A : 



3 



Ao and b) A 



5 



Ao. The 



physical domain admissible for motion is shaded. In the second case there is no physically allowed region for which a > 0. This 
case is distinguished by behaviour at infinity when Hubble function is finite like for the de Sitter solution (see formula (|90|) '). 
This state is the global attractor in the future. All solutions are of the bouncing type. 



where we used definition of the Ashtekar variable A 
(IA3p . From such a particular holonomies we can con- 
struct holonomy along the closed curve a = Dij. This 
curve is schematically presented on the diagram below. 
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This holonomy can be written as 



exp 



^^^V^'^AiA^ere][TurJ,+0{^l') 



(A15) 



where we have used Baker-Campbell-Hausdorff formula 
and the fact that for flat FRW, field strength (|XTO)) sim- 
plifies to the form F^j, = e^^A^A^. Now, equation (jAlSp 
can be simply rewritten to the form 



choose that area Hi to correspond to physical area a'^/i^ 
or to the area /i^. So in the case = V^^^a?^'^ = 
we have the limit 



(A18) 



This approach we call /x— scheme. In the case Hi — |po|m^ 
where po correspond to the eigenvalue 



(A19) 



tr 



-2 lim ■ 



21^2/3 



(A16) 



The trace in this equation can be calculated with use of 
definition (|A14|) 



tr 



Tk 



^kij 



sin^ (lie) . 



(A17) 



In Loop Quantum Gravity the limit fi —^ in the formula 
(|A16|) does not exist because of existence of the area gap. 
The area gap corresponds to the minimal quantum of 
area A = 2\/3Trjlpi, which arises as the first non-zero 
eigenvalue of the are operator [2^ . So instead of the limit 
in equation (|A16P we should stop shrinking the loop at 
the appropriate minimal area A. This area corresponds 
to the area intersected by the loop. For the holonomy 
in the direction °efda the area is as explained in the 
diagram below (where i — k A j). 




The limit ^ A corresponds to fi ^ fi. Now, we 
must connect the area Hi with the length /i. We can 



and taking the limit 



fJ-^ fJ' = fJ-o 



3V3 



(A20) 



which we call /ig— scheme. 

In the quantum version, we can combine equations 
(|A16|) and (|A17p . and write 



^ab - 



JJ^ekij t^a ^b' 



(A21) 







whereas in the classical case we have 

^ab — 2/Z^^^3 a b' 
^0 



(A22) 



So, from these two equations we see that quantum effects 
can be introduced by a replacement 



sin (fic) 



(A23) 



in the classical expressions. The Hamiltonian (|A11[) with 
holonomy correction takes the form 



sin (/ic) 



(A24) 
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